Flat bands with non-trivial topology in three dimensions 
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We construct a simple model for electrons in a three-dimensional crystal where a combination 
of short-range hopping and spin-orbit coupling results in nearly fiat bands characterized by a non- 
trivial Z2 topological index. The flat band is separated from other bands by a bandgap A much 
larger than the bandwidth W . When the flat band is partially filled we show that the system remains 
non-magnetic for a significant range of repulsive interactions, fn this regime we conjecture that the 
true many-body ground state may become a three-dimensional fractional topological insulator. 



Introduction. - When single-particle states of bosons 
or fermions comprise a flat band, the effect of interac- 
tions is generically non-perturbative and can lead to new 
states of quantum matter. For 2D electrons in a strong 
perpendicular magnetic field the flatness of Landau lev- 
els combined with non-zero Chern numbers lead s to the 
spectacular fractional quantum Hall (FQH) effects^! and 
other exotic states such as the Wigner crystalP Recently, 
models for spinless fermions moving in two-dimensional 
lattices have been constructed that exhibit nearly-flat 
bands with non-zero Chern numbers, but remarkably in 
zero external magnetic field.^^ It has been proposed, 
on the bas is of heuristic arguments^' and numerical 
simulations j 6 * 8 * 9 ! that at partial filling repulsive interac- 
tions may drive these systems into the zero-field FQH 
states. Although it is not clear apriori how to experimen- 
tally realize the above situation in a physical system, the 
prospect of engineering systems that can support frac- 
tional excitations without involving a magnetic field has 
led to significant interest recently in flat or nearly flat 
bands in various lattice models. 

In this Communication we ask whether it is possi- 
ble to achieve a flat band with non-trivial topology in 
a system of electrons in three spatial dimensions. It is 
well known that the concept of the first Chern number 
(and the FQH physics) does not generalize to three di- 
mensions. However, recent advances in the theory of 
time-reversal (7~) invariant band insulators have estab- 
lished a uni quely three-dimensional ^-valued topologi- 
cal invarianlPSHUl called vq which, as we argue below, can 
play a similar role. Taking a specific example of electrons 
moving in the 3D edge-centered cubic (perovskite) lat- 
tice we show that a suitable combination of short-range 
hoppings and spin-orbit coupling (SOC) terms can give 
rise to a nearly flat band with a non-trivial Z 2 index 
vq — 1. With the chemical potential inside the bandgap 
such crystal would be a strong topological insulator with 
an odd number of topologically protected gapless states 
associated with all of its surfaces. At partial filling and 
in the presence of repulsive interactions we argue that 
the system may become a three-dimensional 'fractional' 
topological insulator 

Many lattice models are known to support completely 
flat bands. In 2D these include p-orbitals in the hon- 
eycomb lattice}^ as well as s-orbitals in the kagome, 
dice and Lieb latticesP^HSD 313 examples include the py- 



rochlore and the perovskite lattice. ^ * 20 * 21 ! In their sim- 
plest form however the above tight-binding models do not 
constitute a platform suitable for the study of exotic cor- 
related phases for at least two reasons. First, in all cases 
the flat bands touch other dispersing bands, usually in 
a quadratic fashion. The interactions are thus bound to 
mix states from adjacent bands and this is clearly detri- 
mental to the emergence of correlated phases. Second, 
even when the Hamiltonian is modified to open a gap, e.g. 
by adding SOC as in Ref.[20l the resulting flat (or nearly- 
flat) bands are in all known cases topologically trivial. In 
such topologically trivial flat bands the eigenstates can 
be chosen as exponentially localized around lattice sites 
and, at partial filling, interactions will typically select a 
crystalline ground state with broken translational sym- 
metry and not an exotic featureless liquid that underlies 
the FQH effect. As already argued in Refs. [3ff6)20|a nat- 
ural place to seek such exotic phases is in a flat band that 
is topologically non-trivial. 

It is well known that in 2D a non-zero Chern number 
represents a topological obstruction to the formation of 
exponentially localized Wannier states.^ It is precisely 
this obstruction that tilts the balance in favor of FQH 
states in partially filled Landau levels (although Wigner 
crystal phases are still known to arise at very low fill- 
ing fractions). For Z 2 -odd phases of 7~-invariant band 
insulators in 2D and 3D a similar topological obstruc- 
tion exists^ if one insists on Wannier states that respect 
T. This consideration suggests that interacting electrons 
partially filling a flat band with a non-trivial Z 2 index in 
a 3D crystal could either (i) spontaneously break T and 
form a Wigner crystal, or (ii) remain 7~-invariant and 
form one of the proposed 3D fractional topological in- 
sulators characterized either by spin-charge separation^ 
or by fractional magneto-electric effect controlled by frac- 
tional values of the axi on par ameter 6 accompanied by 
ground state degeneracyp3H2D 

In the following we lay groundwork for future investiga- 
tions of these 3D exotic phases by constructing a simple 
tight-binding model whose spectrum has a flat band char- 
acterized by vq = 1 and is separated from other bands 
by a large gap. We remark that constructing a flat band 
with non-trivial topology in a 3D crystal is mathemati- 
cally more constrained problem than in 2D and it is by no 
means obvious that this can be achieved with short-range 
hoppings only. 
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The model. - As a first step in the program outlined 
above we study a simple tight-binding model with intrin- 
sic SOC terms on the 2D Lieb lattice, seen in Fig.[lji, and 
the 3D peroviskite lattice seen in Fig. [2^,. The relevant 
Hamiltonian reads 
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where c ia creates an electron of spin a on site i of the 
lattice, Uj are the hopping amplitudes which we take 
equal to tt, £2, and £3 for the first, second and third near- 
est neighbor sites respectively, and S is an onsite energy 
for the corner sites (indicated as filled circles in Figs, la 
and 2a). Lastly, A represents the amplitude for the next- 
nearest neighbour SOC where d} 4 and d? are the two 
unit vectors along the nearest neighbour bonds connect- 
ing site i to its next-nearest neighbour j and er is the 
vector of Pauli spin matrices. 

In Ref. [201 we showed that for nearest neighbor hop- 
ping (ti only) and non-zero A both Lieb and perovskite 
lattices became topological insulators possessing a non- 
trivial Z2 invariant. In both cases the bands carrying 
non-trivial invariants were strongly dispersing. From 
here then, our goal is to flatten out these bands by tun- 
ing the hopping strengths £2, ^3, A and the onsite energy 
5 without leaving the topological phase. Also, we want 
the flat band to be separated from all other bands by a 
large gap A so that interactions do not mix states from 
different bands. The relevant figure of merit then, is 
the ratio of the flat-band width W to the bandgap A. 

Lieb lattice. - As a warm-up exercise we first consider 
the valence band in the 2D Lieb lattice. Going over to 
momentum space and diagonalizing the Hamiltonian "Hk 
that follows from Eq. , the result of the tuning proce- 
dure can be seen in Fig. [TJd. 

The procedure used for the tuning was the simplest 
available, namely introducing the parameters one by one 
and modifying the values by hand until the band was as 
flat as we could achieve. For the following set of val- 
ues: t x = 1, t 2 = -0.17, f 3 = -0.114, A = 0.303 and 
5 = —0.725, the ratio of the bandwidth to the bandgap, 
A/W « 40.45. A systematic numerical minimization 
technique could probably improve upon our result, but 
comparing to previous results on the kagome, honey- 
comb, square, checkerboard and ruby latticesj^l this 
value is already quite respectable. 

To show that the system has not migrated away from 
the topological phase during the parameter tuning stage, 
we first solve the system in a strip geometry numerically 
using the same parameters as above. This verifies that 
the spin-filtered gapless edge states, seen in Fig. [TJ;, in- 
deed persist. We also consider the Hamiltonian 



H' (x) = (l-x)H + xH, 



(2) 




where Hq includes only t\ and A and is known to support 



FIG. 1: (Color online) (a) Lieb lattice showing 3-site ba- 
sis in unit cell and basis vectors ai and &2- Second and third 
neighbor hoppings are indicated by red and green dotted lines, 
respectively, (b) Tight-binding dispersion for after tun- 
ing parameters to achieve maximal figure of merit A/W. For 
clarity the bulk energy bands are plotted over the shifted Bril- 
louin zone < k x < ir/a, < k v < n/a where a represents 
the length of the nearest neighbor bond, (c) Bandstructure 
for strip of N y — 10 unit cells with open boundary conditions 
along the y direction and infinite along x with the same pa- 
rameters as in (b). (d) Semi-Log plot showing the evolution of 
the gap as one moves between a known topological insulating 
phase and the final set of finely tuned parameters. 



the topological phaseP^The result for A/W, starting at 
A = 0.1 (ti = 1) then adiabatically tuning x from to 
1, for the same final set of parameters in H, can be seen 
111 Fig. [ijd). The ratio increases from its initial value 
A/W « 0.17 to its largest value 40.45 at x = 1, without 
closing the gap. The system is seen to remain in the 
topological phase. 

What can one say about the fate of the ground state in 
the presence of repulsive interactions when the flat band 
is partially filled? Close to half filling, since the spin 
up and down electrons are decoupled, one expects the 
Stoner instability towards the ferromagnetic state. If we 
parametrize the interaction by the usual on-site Hubbard 
U then a gap ~ U should open between the majority and 
minority bands when W < U < A. For large enough U 
one can focus on the spin-polarized electrons in the ma- 
jority band. In the presence of sufficiently strong residual 
interactions (e.g. nearest-neighbor repulsion V) and at 
fractional filling, this problem becomes similar to that 
considered in Refs. [1HS] with spinless fermions. Anal- 
ogous arguments then suggest the emergence of FQH 
states under favorable conditions. Another possibility 
is the formation of a T-invariant fractional topological 
insulatop22 which can be pictured as two decoupled FQH 
states for the two spin species. Such a state might be 
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FIG. 2: (Color online) (a) Perovskite lattice showing 4 sites in 
unit cell along with basis vectors. Second and third neighbor 
hoppings are indicated by red and green dotted lines, respec- 
tively, (b) High symmetry points in the Brillouin zone, (c) 
Bandstructure inside bulk along path of high symmetry for 
"Hit with finely tuned parameter set. (d) Semi-Log plot show- 
ing the evolution of the gap as one moves between a known 
topological insulating phase and the final set of finely tuned 
parameters. 



favored if U < V. 

Perovskite Lattice. - We now approach the main ob- 
jective of this work: constructing flat bands with a non- 
trivial Z 2 index in a 3D lattice. Our starting point is 
again the Hamiltonian given in Eq. Q but now on the 
perovskite lattice, shown in Fig. [5Ja). In Ref. I2TJ1 we 
showed that for nearest-neighbor hopping and A > the 
system becomes a (1; 111) strong topological insulator 
when the lowest doubly degenerate band is filled. In an 
analogous fashion to the two dimensional case then, we 
aim to flatten out the lowest energy band by means of 
including additional parameters t 2 , t 3 and 5. The band- 
structure for the set of parameters that maximize the 
bandgap to bandwidth ratio is shown in Fig. [2b, over a 
path of high symmetry in the Brillouin zone illustrated in 
Fig. (2}j. The following set of values: t x = 1, i 2 = -0.416, 
t 3 = -0.047, A = -0.331 and S = 1.318, yielded the ratio 
A/W fa 17.56. We remark that the ratio was calculated 
over the entire Brillouin zone and in Fig. we have 
simply chosen one possible path here to illustrate the 
result. Including additional parameters in H, such as 
longer range hoppings, could no doubt further improve 
the above figure of merit but we do not pursue this here. 

To ensure that the system remains in the topological 
phase, we adiabatically tune the parameter x from to 1 
in the Hamiltonian H'(x) defined in Eq. The result 
of this procedure can be seen in Fig. [2[i. The ratio, again, 



remains finite across the entire range from its initial value 
A/W pa 0.15 to its largest value 17.56 at x = 1. We 
conclude that the flat-band carries Z 2 index (1;111). 

Now imagine that the (doubly degenerate) flat band 
is partially filled and consider the effect of repulsive in- 
teractions. Unlike the 2D case discussed above where 
the residual U(l) spin symmetry is preserved despite the 
presence of SOC, in the 3D strong topological insula- 
tor the SU(2) spin symmetry is completely broken. In 
this situation it is not clear what the leading instabil- 
ity might be. In a similar flat-band setting describing 
Ir-based pyrochlore Y2&2O7 Pesin and Balents 1 - argued 
for an exotic T-invariant spin-charge separated topologi- 
cal Mott insulator while others found more conventional 
magnetic phasesP^SU Below, we investigate the magnetic 
instabilities of our model in the simplest case with on-site 
repulsion and at exact half filling of the flat bands. We 
show that the system remains non-magnetic when inter- 
action strength is below the critical value which is large 
compared to the bandwidth W. In this regime, there- 
fore, the true many-body ground state could become a 
T-invariant fractional topological insulator. 

Magnetic Instabilities. - We extend the Hamiltonian 
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FIG. 3: (Color online) (a) Phase diagram for Lieb lattice 
(W w 0.034), and (b) perovskite lattice (W w 0.072). Critical 
values U c , Uc for the onset of magnetic order reflect the charge 
distribution on the basis sites. The latter is approximately 
uniform for the Lieb lattice implying U C ~U' C . On perovskite 
lattice most of the charge density is on site 1 causing U c <C U' c . 
(c) Band structure for Lieb lattice with U/W = U' /W = 15, 



,(3) 



m v ~' ~ 0.35 and m' 2 ' = m 
lattice with U/W = 5, U'/W = 2x n, 
m (3) w 0.11 and m (4) w 0.05. 



0.32, and (d) perovskite 
« cs 0.62, = 
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(1) above by including a Hubbard term 

H^ = Uj2r$r$ + U>j:4!n$ (3) 

i i,l>l 

with the goal of mapping out the magnetic phases as 
a function of U and U' by means of a standard mean- 
field calculation. The superscript I denotes the basis sites 
in the lattice, Figs, la and 2a. We use the following 
decoupling 

riifnn -t n^{n i]r ) + miimf) - (n it ) (n 4 ) (4) 

and definine the magnetization rap = ) — (nf) ) and 

occupation number fif = (nfJ)-\-(nf)). Furthermore, we 

focus on uniform magnetic phases, such that mfp = 
independent of the site index i. 

To map out the phase diagram, we minimize the 
ground state energy with respect to mS 1 '. The result 
for half filling of the flat bands can be seen in Figures [3^ 
and|3j3, for the Lieb and perovskite lattices respectively. 
Here m = mP"' is the total magnetization per unit 
cell. The key observation to be made here is that both 
systems remain non-magnetic over a significant range of 
interaction strengths. 

Conclusions. - We have established that it is possible 
to engineer nearly flat bands characterized by non-trivial 
Z2 topological invariants in both the 2D Lieb lattice and 
its three dimensional counterpart, the edge centered cu- 
bic (perovskite) lattice. With our 2D example we have 
simply added to the growing number of lattices capa- 
ble of producing this behavior, whereas the 3D result is 
to the best of our knowledge completely novel. What 
makes these findings potentially broadly relevant is that 
they provide a concrete framework for future studies ad- 
dressing the nature of the many-body ground state that 
occurs in the presence of repulsive interactions and at 
partial filling. 



Unlike the 2D case where the Laughlin liquicP fur- 
nishes a well established paradigm for the topologically 
ordered correlated ground state, in 3D such a template 
is presently lacking. Therefore, understanding the fate 
of the ground state of electrons in the topologically non- 
trivial 3D flat band in the presence of strong interactions 
beyond the simple mean-field analysis is clearly a prob- 
lem outside the scope of this study. In this work we made 
a first step in this direction by demonstrating, within the 
mean-field approximation, that the ground state remains 
non-magnetic over a significant portion of the U-U' phase 
diagram (Fig 3b). Since mean-field calculations tend to 
overestimate the stability of ordered phases the actual 
non-magnetic region will likely be somewhat larger and 
can in principle host a T-invariant fractional topologi- 
cal insulator. Our result in 3D thus provides a concrete 
foundation for addressing these interesting questions us- 
ing the suite of analytical approaches or numerical many- 
body techniques and we hope that it might inspire further 
studies. 

Finally we remark that there exist many perovskites in 
nature. Our theoretical results demonstrate that a near- 
ideal situation (in terms of strong interaction physics) 
can arise in a very simple model in this family of lat- 
tices. Also alluded to previouslj^ was the possibility of 
engineering a relevant 2D system by modulating a two- 
dimensional electron gas with a periodic potential having 
Lieb lattice symmetry, in an analogous fashion to the 'ar- 
tificial graphene' created recently.^ Another possibility 
lies with cold fermionic or bosonic atoms in optical lat- 
tices as discussed in Refs. |28l29j . with a detailed study 
on the Lieb lattice provided in Ref. [50] . 
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